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Abstract 

We study the complex-temperature properties of a rare example of a statistical mechanical 
model which is exactly solvable in an external symmetry-breaking field, namely, the Ising 
model on the square lattice with jSH = ±m/2. This model was solved by Lee and Yang 
|]1[]. We first determine the complex-temperature phases and their boundaries. From a 
low-temperature, high-field series expansion of the partition function, we extract the low- 
temperature series for the susceptibility x to 0(u^^), where u = e~^^ . Analysing this series, 
we conclude that x has divergent singularities (i) at m = = — (3 — 2^/^) with exponent 
7g = 5/4, (ii) at u = 1, with exponent = 5/2, and (iii) at n = = —1, with exponent 
7^ = 1. We also extract a shorter series for the staggered susceptibility and investigate its 
singularities. Using the exact result of Lee and Yang for the free energy, we calculate the 
specific heat and determine its complex-temperature singularities. We also carry this out for 
the uniform and staggered magnetisation. 



* email : vmatveev@max. physics . sunysb . edu 
* * email: shrock@max .physics . sunysb . edu 



1 Introduction 



The Ising model has long served as a prototype of a statistical mechanical system which 
undergoes a phase transition with associated spontaneous symmetry breaking and long range 
order. In the absence of an external magnetic field H, the free energy of the d = 2 (spin 
1/2) Ising model was first calculated by Onsager 0, and the expression for the spontaneous 
magnetisation first calculated by Yang both for the square lattice. The model has never 
been solved in an arbitrary external magnetic field. However, in one of their classic papers, 
Lee and Yang [|l| did succeed in solving exactly for the free energy and giving an exact 
expression for the magnetisation of the Ising model on the square lattice for a particular 
manifold of values of H depending on the temperature T, given by 

h='im: (1.1) 

2 ^ ^ 

Although this is not a physical set of values, owing to the imaginary value of H and the resul- 
tant non-hermiticity of the Hamiltonian, this model is nevertheless of considerable interest 
as a rare example of a statistical mechanical model for which one has an exact solution in 
the presence of a symmetry-breaking field. Further work on the derivation of the Lee- Yang 
solution was reported in Refs. |^, |^, |^, |^. 

In the present paper, we shall investigate this model in a wider context, generalising the 
temperature to complex values. There are several reasons for studying the properties of 
statistical mechanical systems with the temperature variable generalised to take on complex 
values. First, one can understand more deeply the behaviour of various thermodynamic 
quantities by seeing how they behave as analytic functions of complex temperature. Second, 
one can see how the physical phases of a given model generalise to regions in appropriate 
complex-temperature variables. Third, a knowledge of the complex-temperature singularities 
of quantities which have not been calculated exactly helps in the search for exact, closed- 
form expressions for these quantities. This apphes, in particular, to the susceptibility of the 
present model, which, like that of the zero-field Ising model, has never been calculated. 



2 Generalities and Complex- Temperature Phases 

In this section we shall work out the complex-temperature phases and their boundaries. We 
begin with some definitions and notation. Recall that the Ising model is defined by the 
partition function Z = J2{an} with the Hamiltonian 

n=-J J2 ^nan'-Hj2ar, (2.1) 

<nn'> n 
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where cJn = ±1 are the spin variables on each site n of the lattice, (5 = {ksT)'^, J is the 
exchange constant, < nn' > denote nearest-neighbor pairs, and the units are defined such 
that the magnetic moment which would multiply the H J2n '^n is unity. We shall concentrate 
here on the square (sq) lattice. We use the standard notation K = pj, h = PH, v = tanhK, 
z = e~^^ , u = = e~^^, w = 1/u, and /i = e~^'^. Note that v and z are related by the 
bilinear conformal transformation 

1 — V , ^ 

z = 2.2 

l + v ^ ' 

It will also be useful to introduce two elliptic moduli. The first is 

which occurs in elliptic integrals in the exact expressions for the internal energy and specific 
heat, where we use the abbreviations 



We record the symmetry 



The second elliptic modulus. 



C = cosh(2K) (2.4) 
S = sinh{2K) (2.5) 

u — i> 1/u =^ ft — > K (2.6) 



k - ' - (2 7) 

< 5(^2 + 2)1/2 (1_^)(1+ 6^ + ^2)1/2 

occurs in a natural way in the magnetisation.^] 

The reduced free energy per site is / = —f3F = limi\f^_^ In Z (where A^^ denotes 
the number of sites on the lattice). In addition to the susceptibility itself, it will also be 
convenient to refer to the reduced susceptibility x = P~^X- 

We begin by discussing the phase boundaries of the model as a function of complex 
temperature, i.e. the locus of points across which the free energy is non-analytic. As noted 
in Ref. there is an infinite periodicity in complex K under the shift K K + nin, 
where n is an integer, and, for lattices with even coordination number q, also the shift 
K — >■ {2n + l)i7r/2, as a consequence of the fact that the spin-spin interaction aiCj in Ti is 
an integer. In particular, there is an infinite repetition of phases as functions of complex 
K] these repeated phases are reduced to a single set by using the variables v, z, or u (or 
variables based on these). 

^Note that these differ from the respective ehiptic moduH kq and fc< g which occur in the internal energy, 
specific heat, and spontaneous magnetisation for the Ising model on the square lattice at h — 0. 
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We also note an elementary symmetry involving h for the (spin 1/2) Ising model on a 
general lattice A. The low-temperature, high-field expansion of Z has the form 



Z = e(''/2)^^^e^^'^Z, (2.8) 

where 

Zr = l + T.^ltl^y'^ (2-9) 

n,m 

where the only property of Zr that we need is the fact that it is a polynomial in z and fi. In 
( p.8|) , we assume periodic boundary conditions, but for the free energy, in the thermodynamic 
limit, this is not essential. Also, parenthetically, we note that for a lattice with even q, only 
even powers of z occur in Z^, i.e. Z^ is a series in u, but we shall not need this fact here. 
Now 

h h + niTC =^ yU — > yU (2.10) 

where n is an integer. Hence, under such a shift, the only change in Z is in the prefactor, 
gAfs/i_ Equivalently, in the corresponding (reduced) free energy 

f = {q/2)K + h+\im N-\l + al^rl^"- D (2.11) 

n,m 

the only change is in the second term, h. Therefore, aside from this term, one may, and we 
shall, restrict to the range 

-Y< Imih) < I (2.12) 

without loss of generality. In the present context, we shall consider just the value h = 171/2; 
our results will apply in the same way to h = —mj2. 

It is useful to review the connection between the square-lattice Ising model with h = m/2 
and the Ising model on the square lattice in zero field . This is done by first considering 
the Ising model on the checkerboard (also called generalised square) lattice, defined by 
assigning different couplings Kj, j = 1,...,4, to the bonds of the square lattice, as shown 
in Fig. 1. Again, for discussions of the partition function, we assume periodic boundary 
conditions. The free energy |^ and spontaneous magnetisation [|TU|, are known for the 
zero-field checkerboard Ising model. Now recall the identity 

e'"' = cosh/i + asinh/i (2.13) 

for a = ±1. For h = i'7i/2, this reduces to e^""" = ia^, and hence exp(/i X^n CTn) = Dnl^^n)- 
Next, consider a dimer site covering of the checkerboard lattice, where by site covering, we 
mean that each site is the member of one (and only one) dimer. As is clear from Fig. 1, a 
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simple covering of this sort is provided by each of the bonds of a single type, say those with 
the couplings. We may thus associate pairs of the cr„'s in the above product with the 
dimers of this covering. To do this, we separate one of the two factors of i for such pair and 
place it in front of Z . One then has 

Zch = l""''^ E ( n {^^rCTs)) exp( E ^nKnn'CTn') (2.14) 

{o-„} <n,n'> 

where ch denotes checkerboard and Knn' refers to the appropriate Kj, j = 1, 2, 3, 4 depending 
on which bond connects the sites n and n' (c.f. Fig. 1). Then one can use the identity ( glSD 
again to write 

Z = z^^/2Eexp( E ^ni^;„'^n') (2.15) 

Wn} <nn'> 

that is, 

Z{{K,}; h = = {tf'^'ZUK'^y, h = 0),, (2.16) 
where K- = Ki^ 2 = 1,2, 3, and 

777- 

K = K, + - (2.17) 

Hence, for the free energy, 

771 

f{{K,}- h = zn/2U = ^ + fHKiy, h = OU (2.18) 

Then, setting Ki = K, i = 1, 2, 3, 4, one can obtain the Lee- Yang result for f{K, h = m/2) 
|l[] from the (analytic continuation of the) free energy for the zero-field checkerboard lattice 
0]. The same method works for the magnetisation and yields the relation 

M({ir,}, h = it: 12),,, = M{{K^}, h = 0),h (2.19) 

and for the m-point correlation functions, which satisfy 

<(Jm--- (Tn^ > {{Ki}, h = nr/2)ch = < (T„i ■ ■ ■ (T„^ > ({K-}, h = 0)ch (2.20) 

Further, it follows from the special case of ( |2.20| ) for 2-spin correlation functions, together 
with the expression for the susceptibility as a sum over the connected (conn.) 2-spin corre- 
lation functions, 

X = E < ^Of^r >conn. (2.21) 
r 

(where < aoaj. >,onn.=< o"oO"r > — M^) that 

X({i^.}, h = 17, /2U = X{{K[}, h = (2.22) 
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Of course, the Ising model with zero field is not equivalent to one with nonzero field, since 
in the former case the partition function and free energy are exactly invariant under the 
transformation cr„ — > — cr„, whereas in the latter case this symmetry is broken explicitly by 
the external field term. This inequivalence is manifested in the fact that eqs. ( p.l6 )-( p.l7 ) 
are not Z2-invariant. Thus under the transformation cr„ — which is equivalent to 
h —h, the i's in these equations are replaced by —i. It is also manifested in the fact that 
while the zero-field Ising model always has a Z2-symmetric paramagnetic (PM) phase, the 
model with nonzero external field h does not have any PM phase. The usefulness of eq. 
( p.l8| ) stems from the special feature that for h = ±z7r/2, this non-invariance is localised to 
just a constant term in the free energy. 



The (reduced) free energy is 



fiK,h = ±in/2) 



±— + n2- 
2 



where C and S were defined in 
the symmetries 

U{P,J,H 



4 J-^ J-n (27r)2 

^ and 

in 
2P 



In 



C^ + 5^-l + S2(cOs(01+02 



-COS 



From U 



) = UiP,J,H 



(2.23) 

—df/d(3 = —Jdf/dK, one has 



ITT 



) 



Ui/3,-J,h=- 



17C 



U{-f3,J,h = -) = -U{P,J,h="-^) 



(2.24) 
(2.25) 
(2.26) 



Similarly, from C = ksK'^d'^f/dK'^, one has 



C{K,h 



and 



C{K,h='-^] 



C{K,h 



2 ^ 



(2.27) 



(2.28) 



The free energy is trivially divergent at = ±cxd, i.e. m = 0, oo; however, this will not be 
important here since these are isolated points and not part of any phase boundaries. The 
curves along which the free energy is non-analytic are given by the locus of points where the 
argument of the logarithm in the integrand of eq. (|2.23| ) vanishes. Expressed in terms of the 
variable u, / is 



?7r 1 

f{K,h = ±m/2) = ±-+-\n 



:i-u) 



1 

+ 4 



deide2 
W 



In 



--2mP(0i,02) 



(2.29) 



^The Hamiltonian in Ref. |i| was defined with a different zero point of tlie energy than that used here. 
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where 

P(^l,^2) = COS^l +008^2 (2.30) 

The above locus of points where the argument of the logarithm vanishes is given by the 
solutions of the equation 

{l + uY-2ux = (2.31) 

where x = P{6i,62), taking values in the range —2 < x <2. These are integrable singular- 
ities. Since the coefficients in this equation are real, the solutions are either real or consist 
of complex conjugate pairs. Moreover, under the replacement u — 1/u, eq. ( p.31|) retains 



its form, up to an overall factor of u^'^. Consequently, the locus of solutions is also invariant 
under this mapping u 1/u. The solutions are shown in Fig. 2(a) and consist of the union 
of the unit circle 

u = e''^ , -TV <(f)<7r (2.32) 



and the finite line segment 
where the inner endpoint is 



e 



1 

— <U<Ue (2.33) 



-(3 - 2^2) = -0.171572875... (2.34) 



and the outer endpoint is l/u^ = —(3 + 2^/2) = —5.828427.... Note that Ue = —Uc, where 
Uc is the usual critical point in the zero-field square lattice Ising model separating the Z2- 
symmetric, paramagnetic (PM) phase from the phase in in the Z2 symmetry is spontaneously 
broken by long-range ferromagnetic (FM) long-range order. 

It is of interest to see how the solutions to eq. (^.31|) are traced out in the complex u 



plane as x varies. For x = 2, this equation has a double root at m = 1. As x decreases 
from 2 to 0, this root splits into a complex conjugate pair, the members of which move 
counterclockwise and clockwise along the unit circle, and finally rejoin to form a double root 
at M = — 1 when x = 0. As x decreases from to —2, this double root again splits, but 
this time into two reciprocal real roots, one of which moves to the right, from u = —1 to 
the endpoint Ue and the other of which moves leftward to m = l/ug. The corresponding 
phase boundaries in the z plane consist of the unit circle = 1 together with the two line 
segments from z = ±Ze = ±i{\/2 — 1) upward and downward along the imaginary axis to 
z = ^1/ Zf. = ±i{\^ + 1), respectively. 

The circle ( ^.32| ) divides the u plane into two separate phases. A fundamental property of 
this model is that the nonzero external field breaks the Z2 symmetry explicitly, so that there 
is no Z2-symmetric phase. Even without using the known expression for the magnetisation. 
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one can identify the phases in this diagram as follows. For sufficiently large real K, the 
interaction of the external magnetic field with the spins is negligible compared with the 
spin-spin interaction, which thus produces a ferromagnetically ordered phase, just as it does 
in the model with /i = 0. This shows that the neighborhood of the origin in the u (or z) plane 
is ferromagnetically ordered. By analytic continuation, it then follows that the entire region 
inside the unit circle |m| = 1 is a ferromagnetically ordered phase, and this is so denoted in 
Fig. 2(a). Similarly, for sufficiently large negative K, the interaction of the external field 
with the spins is again negligible compared with the spin-spin interaction, which produces a 
phase with antiferromagnetic (AFM) long-range order. By analogous analytic continuation 
arguments, it follows that the entire region outside the unit circle is the AFM phase. This 
may be shown as follows. One may think of the phase diagram in the complex variable 
w = 1/u. By the u ^ 1/u symmetry noted above, the phase boundaries in this variable 
are the same as those in Fig. 2(a). The above argument shows that the neighborhood of 
the origin is antiferromagnetically ordered, and hence, by the same analytic continuation 
method, the entire region inside the unit circle |w| = 1 is an AFM phase. Mapping this back 
to the u plane, one has shown that the entire region outside of the unit circle |u| = 1 is the 
AFM phase. Since the u plane consists of precisely two regions with the property that within 
each one can analytically continue from any one point to any other, we have thus obtained 
a complete description of the phases in the model. The same reasoning implies that in the z 
plane, the phase with \z\ < 1 is FM and the phase with \z\ > 1 is AFM. In passing, we note 
that complex-temperature properties of the h = Ising model on ci = 2 lattices have been 
studied in Refs. 0, 0, |14|, |15|, ||, 0, |18|, [19|, |2^. In the case /i = for the square lattice. 



the analogous locus of points across which the free energy is singular form a limagon ||18i 
defined by Re{u) = 1 + 2^/^ cos + 2 cos 2uj, Im{u) = 2^/^ sin + 2 sin 2cj for < < 27r, or 
equivalently, in the z plane, the circles [jl2| z = ±1 + -\/2e*^, for < 6' < 27r. 

By the use of the conformal mapping (|2.2|) on the z plane, or by re-expressing the free 
energy in terms of the variable v and again solving for the locus of points where the argument 
of the logarithm vanishes, we find that the phase diagram of the model in the v plane is 
as shown in Fig. 2(b). The unit circle \z\ = 1 is mapped to the imaginary axis in the v 
plane, and the respective line segments from z = ±Ze to are mapped to the arcs from 

V = e^™l^ to f = e^'^'^'"!'^. It is interesting that the Re{v) = (i.e. imaginary) axis forms the 
boundary between the complex-temperature FM and AFM phases. One may understand 
this by recalling that (i) if h were real and positive (negative), this would favor FM (AFM) 
ordering, but a pure imaginary value of h does not favor FM over AFM ordering, or vice 
versa; (ii) similarly, if the spin-spin coupling K is real and positive (negative), it favors FM 
(AFM) ordering, but a pure imaginary value of K (and hence v) does not favor FM over 
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AFM ordering, or vice versa. Therefore, if both h and K are pure imaginary, as is the case 
here for the imaginary axis in the v plane, then the system is precisely balanced between FM 
and AFM order, so that this axis should be the boundary between the complex-temperature 
FM and AFM phases, and this is just what our explicit calculation shows. 

The mapping defined by m — k^, where k was defined in (|2.3| ), takes the the locus of 
points ( p.32|) and (|2.33| ) to a single semi-infinite line segment extending from 1 to oo in the 
complex plane. All points in the plane are analytically connected to all other points. In 
particular, the mapping u ^ k"^ takes both the complex-temperature FM and AFM phases 
in the u plane to the same respective regions in the plane, as is clear from the symmetry 
( p.6|) and the fact that the transformation u ^ 1/u interchanges the FM and AFM phases 
in the u plane. 



3 Complex- Temperature Behaviour of the Internal En- 
ergy and Specific Heat 

3.1 Exact Expressions 



From the free energy ( |2.29| ), it is straightforward to calculate the internal energy U and 
specific heat C (per site). In terms of the variable u, we find that 



U = -J 



l + u 



u 



u 



l + u 



(3.1.1) 



where the elliptic modulus k was given above in eq. (|23) and i^(A;) = /o"/^(l-P sin G) ' d6 
is the complete elliptic integral of the first kind. This expression holds for both the FM 
and AFM phases and exhibits the symmetries ( |2.25| ) and ( |2.26| ). Since either of these has 
the effect of taking u 1/u, and since this mapping takes the interior of the complex- 
temperature FM phase to the interior of the complex-temperature AFM phase, the values 
of U in these two phases are simply related by ( |2.25D -( p.26D . In the FM phase, the first few 
terms of the small- |m| expansion (complex-temperature generalisation of the low-temperature 
expansion) are 



U 



-2J 



1 + - I2u^ + 60m^ - 280m^ + 0{u^ 



(3.1.2) 



In the AFM phase, the corresponding expansion parameter is w = 1/m, and U has the same 
expansion as ( |3.1.2| ) with J replaced by — J and u replaced by w. 

As discussed above, in the limit J oo, and hence i^' ^ oo with fixed H , the spin-spin 
interaction overwhelms the contribution of the external field coupling, which therefore has 
a negligible effect, to leading order. It follows that in this limit, the value of the internal 
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energy should be the same as the value for h = 0, i.e., 

U{u = 0,h = m/2) = U{u = 0,h = 0) 



(3.1.3) 



It is interesting to compare the small- |m| series expansions of these two functions to ascertain 
the finite-u corrections to this equality. For this purpose, we recall that B 



U{K,h = 0) = -J 



1-u 



where 



4^(1 -m) 



(3.1.4) 



(3.1.5) 



{1 + uy 

The expression (|3.1.1CI| ) holds for all phases, PM, FM, and AFM. In the FM phase, it has 
the small- |m I expansion 



U{h = 0) 



-2J 



(3.1.6) 



Clearly, the expansions ( |3.1.2| ) and ( |3.1.6| ) agree with the relation ( ^.1.3] ) for u = 0. U {K, h = 
0) also satisfies the symmetries analogous to ( |2.25| ) and ( |2.26D , with h = z7r/2 replaced by 
/i = 0; as a consequence, in the AFM phase, the small-|w| expansion of U{K, h = 0) is given 
by ( p.l.6| ) with J replaced by — J and u replaced by w. 

For C we get 

^2 (1 + ^2^ 



c 



u 



u 



(1-m)2 7r(l + 6M + M2) 



7r(l + M) 



(3.1.7) 



where K{k) was defined above and E{k) = Jq^'^{1 — k'^ srn^ OY^'^dO is the complete elliptic 
integral of the second kind. Again, this expression holds for both the FM and AFM phases. 
It will also be useful to express C in an equivalent form, using ( p.3|) : 

{1 + u''){1 + k)K{k) 



c 



u 



EiK) 



[l-uf 7r(H-fi;) 



7r(l + Qu + u^) 



(3.1.8) 



C/K"^ has the small- |m| expansion 



C 



-64m^ + 288m^ - 1920m^ + 11200m^ + 0{u^ 



(3.1.9) 



^keK"^ 

For comparison, the specific heat for the Ising model on the square lattice with /i = is 
C 4(1 -k') 



2{K{ko) - E{ko)} - (1 - + k',K{ko)} 



(3.1.10) 
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which has the small- !«! expansion 

C 



6Au' + 288^^ + 1152m^ + 4800m^ + 0(m'') (3.1.11) 



Of course, in this case, the positivity of the specific heat requires that the coefficient of the 
lowest order term must be positive (the first negative coefficient occurs in the term). We 
proceed to determine the complex-temperature singularities of U and C for the present case, 
h = m/2. 

3.2 Vicinity of w = 

As discussed in connection with Fig. 2(a), the point m = is the endpoint of the singular line 
segment protruding into the complex-temperature extension of the FM phase. All approaches 
to this point, except directly from the left along this singular line segment, occur from within 
the complex-temperature FM phase. As m — »• Me = 

-(3-2^/2), K ^ -1. The internal energy 

diverges like 

U ^/2 

— ^ — ln(l — w/ue), as u ^ Ue (3.2.1) 

J 71 

In the specific heat, the dominant divergence arises from the term in ( |3.1.7| ) involving E{k,) 
and is 

fc^^- vr(l-nK) ' (3.2.2) 
so that the associated singular exponent for C at m = is 

<,FM = = 1 (3.2.3) 

where the subscript FM indicates the phase from which this point is approached, and the 
prime is the standard notation indicating that the approach to this singular point is from 
within a broken-symmetry phase. (The last feature is, of course, true of all of the singular 
points for h ^ 0.) There is also a weaker, logarithmic divergence arising from the term 
involving K{k). The value of K at Ue in ( |3.2.2| ) is 



Ke = —- InUp 
4 



1 , 1 



ln(3 - 2^/2^ +m + 2nin] (3.2.4) 



where n labels the Riemann sheet used for the evaluation of the logarithm, which we shall 
take to be n = below, unless otherwise indicated. 
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3.3 Vicinity of u = 1/ue 

The point m = 1/ue is the left end of the singular line segment protruding into the complex- 
temperature extension of the AFM phase. Except for the approach directly from the right 
along the singular line segment, all approaches to this point occur from within the complex- 
temperature AFM phase. As m — > l/wg, /t — > —1, as is clear from the previous remarks and 
the symmetry (|2.6| ). The internal energy again diverges like 

^ ^ 1 (',',^^ 
> m(l — Uf,u) , as u ^ — (3.3.1) 

J TT Ue 

In the specific heat, the dominant divergence again arises from the term in ( ^.1.71 ) involving 
E{k) and is 

C 4V2 

^ ^ nil-ueu) ' " ^ '/"^ 

so that the associated singular exponent for C at the outer endpoint (oe) u = 1/ue is 

O^'oeAFM = Oi'oe=^ (3.3.3) 



The value of K corresponding to m = l/ug in eq. ( 3.3.2| ) is, for the principal Riemann sheet 
of the log, Koe = -(l/4)[ln(3 + 2^/2) + in]. 

3.4 Vicinity of w = — 1 

As M ^ — 1 (denoted Us), k diverges; if we set u = —1 + ee^'^ and let e — > 0, then k ~ 
— 4e~^e~^*'^. One easily sees that U{u = —1) is finite. For C, we observe that the first and 
second terms on the right-hand side of eq. ( |3.1.7| ) or ( |3.1.8| ) are finite. By the use of the 



elliptic integral identity (1 + k)K(k) = K{2k^/'^ / [1 + k)) (see, e.g. [^), we can rewrite the 
term involving K{k) in eq. ( p. 1.8] ) as 

<l+"^' -K(—) (3.4.1) 



As M — > — 1 and k diverges, K{2k^/'^ / [1 + k)) -^(0) = 7r/2, so that C is finite, although 
non-analytic, at m = —1. This is true for the approach to m = — 1 from either the FM or 
AFM phases. We thus have 

«s,FM = «Ufm = {log. finite) (3.4.2) 
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3.5 Vicinity of w = 1 



Asm— ^-IjK^l. In U the leading potential singularity arises from the first term in ( |3.1.1| ) 

-2 



U 
1 



as u 



u 



(3.5.1) 



Now K = — (l/4)lnM, so that, if one uses the first Riemann sheet of the logarithm, then 
M — > 1 maps to ^ 0. Recalling that K = jSJ, if this zero in K is due to /3 — at fixed 
nonzero J, then eq. ( p.5.1| ) shows that U diverges for u ^ 1; however, if the zero in K is 
due to J —> at fixed nonzero /5, then, expanding ( p.5.1| ), one finds that U — > 1/(2/3). 
For the specific heat, from ( p.l.7| ), it follows that as m — 1, 



(i-uy 



1, / 32 
H n 

477 Ml 



- u 



+ 



(3.5.2) 



where ... refers to less singular terms. If we again use the principal Riemann sheet of the 
logarithm, so that u ^ 1 corresponds to i^' ^ 0, then (|3.5.2|) becomes 



as u 



1 



(3.5.3) 



2 ■ TT ^K^' 

That is, C has a finite logarithmic singularity at this point, and hence a corresponding 
exponent 



a 



l,FM 



a 



1,AFM 



{log. finite) 



(3.5.4) 



In passing, we note that if one were to use a Riemann sheet different than the principal 
(n = 0) one in evaluating K = —(1/4) ln(l), so that K = —inn/2 ^ 0, then C would diverge 
quadratically at u = 1. 



3.6 Elsewhere Along the Singular Curves 

We discuss here the behaviour of U and C as one crosses the singular locus of points comprised 
by the unit circle ( p.32| ) and the line segment ( p.33|) away from the points u = Ue, 1/ue, — 1, 
and 1. The singularities which one encounters in this case are associated with passage across 
the branch cut of the elliptic integrals in (|3.1.1| ) and (|3.1.7|) . We recall that the elliptic 
integrals K{k,) and E{k) are analytic functions of in the complex plane except for 
respectively divergent and finite branch points at = 1 and an associated branch cut, 
which is normally taken to run from k'^ = 1 to k"^ = oo along the positive real axis in this 
plane. To illustrate the nature of the singularities, we shall consider moving outward along a 
ray in the u plane defined hj u = pe*^ with p increasing from to oo at fixed 6, say 6 = n/G. 
As shown in Fig. 3, the image point in the plane also moves out from the origin, starting 
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with an angle of vr/S but bending around to the right. As we cross the unit circle in the u 
plane, leaving the FM phase and entering the AFM phase, the image point in the plane 
crosses the branch cut moving vertically downward. This branch cut is precisely the image 
of the unit circle |m| = 1 (and also of the singular line segment (p.33 ).) For 6 = 7r/6, the 



crossing point is at = 2'^{7 — 4\/3) = 1.1487.... In the plane, one thus passes onto the 
second Riemann sheet of the elliptic functions K{k,) and E{k). If one projects back to the 
first Riemann sheet, these functions have discontinuous imaginary parts across this branch 
cut. As p continues to increase toward oo, ~ 16p~^e~^^ so that the image point curves 
around finally approaches the origin in a "northwest" direction, at an angle of — 7r/3, but on 
the second Riemann sheet. In Fig. 3 we show the image point for p in the range from to 
30. 

4 Complex- Temperature Behaviour of the Uniform and 
Staggered Magnetisation 

The magnetisation M is ||T], Q 

Miu,h = m/2) = ' -— (4.1) 

Note that {l + Qu + u^) = {l-u/ue){l — Ueu). By analytic continuation, this formula holds 
throughout the complex-temperature extension of the FM phase. The identity discussed 
above, and the resultant eq. ( p.l9| ) yields the relation 



M(m, h = ±tn/2) = M{-u, h = 0)"^ (4.2) 

where 

, , , (1 + U)1/4(1 -6^ + ^2)1/8 

M{u,h = 0)=^- ^-f-;^, (4.3) 

As is well known, one can express M{u, /i = 0) as 

M(«,/i = 0) = (l-(fc<,o)^)V« (4.4) 

where ^ ^ 

= sinh2(2ir) = (T^ ^^-^^ 
This quantity also enters in exact expressions for correlation functions in the FM phase of 
the h = Ising model [p2| , p^ , [2^ . Given the relation ( [4. 21 ), it is natural to write 



M(M,/i = m/2) = (l-(A;<)2)i/8 (4.6) 
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where the elhptic modulus was introduced in eq. ( p.7|) . The magnetisation for the h = i7i/2 
case vanishes continuously at the point u = —1 (denoted Ug) with exponent 



diverges aX u = with exponent 
and diverges at m = 1 with exponent 



Ps = \ (4.7) 

/?e = -\ (4.8) 

A = -\ (4.9) 

Elsewhere on the boundary of the complex-temperature extension of the FM phase, i.e., the 
unit circle in the u plane, M vanishes discontinously. Note that the apparent divergence at 
the point u = 1/ue does not actually occur, since this is outside of the complex-temperature 
FM phase, where the above analytic continuation is valid. 

The staggered magnetisation Mgt does not seem to have been explicitly discussed in the 
literature, but one can easily obtain it, as follows. Mat may be defined via 

Ml^ = lim < aoar > (4.10) 

r|-— >oo 

where 

ar = (-l)^(^Vr (4.11) 

where 

2 

p{r) = Y.ri (4.12) 

1=1 

i.e. 5"r = (Tr for r on the same sublattice as r = and —ay. for r on the other sublattice of 

the (bipartite) square lattice. To evaluate Mst via eq. ( |4.1CI| ), it suffices to take r = (r, 0) 



or (0,r), i.e. the row or column 2-spin correlation function. From the known asymptotic 
behaviour of this correlation function one immediately finds that 

M,t{w) = M{u w) (4.13) 



where, as before, w = 1/u. This is consistent with eq. (f4.6|) since (c.f. (|2.6| )) m — > 1/m takes 
fc< — fc<, and /c< enters squared in ( [4.6|) . Of course, Mgt vanishes identically outside the 
complex-temperature extension of the AFM phase. Further, we may immediately conclude 
that Mst vanishes continuously at u = —1 with exponent ( |4.7|) , diverges at u = l/ue with 
exponent (|4.8|) , and diverges at m = 1 with exponent ( [1.9|) . Elsewhere along the bound- 
ary of the complex-temperature AFM phase, Mgt vanishes discontinuously, with the same 
discontinuity as M. 



14 



It is of interest to compare these results with the behaviour of M and Mst for /i = 
(again on the square lattice). Aside from the physical PM-FM and PM-AFM critical points 
u = Uc = (3 — 2^/^) and 1/uc, where, respectively, M and Mst vanish continuously with 
exponent (3 = 1/8, they also both vanish continuously at the complex-temperature point 
u = —1, with the same exponent, (3 = 1/4. Note that for h = there is only one point, viz., 
u = —1, where the FM and AFM phases are contiguous and M and Mst vanish continuously, 
whereas for h = 111/2 there are two such points, namely, u = —1 and u = 1. 



5 Extraction and Analysis of the Low- Temperature Se- 
ries for X 

5.1 Generalities 

In order to investigate the complex-temperature singularities of the susceptibility x, we 
shall make use of the low-temperature, high-field series expansion for the free energy or 



equivalently the partition function of the Ising model on the square lattice I^S], |26|, In 
Ref. Baxter and Enting calculated this expansion for the partition function to order 

0(m^^). The series for Zr in eq. (|2.9|) is 



n=2 m 



where j < m < j"^ for n = 2j and j < m < j{j — 1) for n = 2j — 1. We extract the series for 
h = ire /2 by calculating x = d'^f/dh'^ and then subsituting fi = —1. This has the form 

00 

xiu,h = t7r/2)=Au\Y.^nU'') (5.1.2) 

n=0 

The results for the c„ are listed in Table 1; the series for Z and the resultant series for x to 
0(m^^) yields the c^'s to order n = 21. Parenthetically, we note that x has been calculated to 
0(m^^) in Ref. and to 0{u^^) in Ref. but the low-temperature, high-field expansion 
of the partition function as a function of yU = e~^'* which would be necessary to extract 
x{h = iTx/2) was not given in these papers (it would be a rather long expression). 

We have analysed this series using dlog Pade and differential approximants. For a recent 
review of these techniques, see Ref. Our notation for these approximants follows 

Ref. and our earlier work on complex-temperature properties of the h = Ising model 
IjT^, |Ty], |2^. In particular, we use first order differential approximants (i.e., i^' = 1 in our 
previous notation); as before, we used unbiased approximants so as to be able to use an 
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n 


Cn 





-1 


1 


8 


2 


-48 


3 


304 


4 


-1863 


5 


11368 


6 


-68840 


7 


414872 


8 


-2490437 


9 


14903648 


10 


-88963696 


11 


529939176 


12 


-3151205475 


13 


18710180192 


14 


-110948037424 


15 


657164715520 


16 


-3888670886593 


17 


22990566432904 


18 


-135819110416784 


19 


801806651588848 


20 


-4730485389238263 


21 


27892958533539784 



Table 1: Low-temperature series expansion coefficients for h = m /2) in eq. (|5.1.2D. 
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extrapolation method for extracting critical exponents. Since the prefactor is analytic, 
we have actually performed the analysis on the reduced function 

X. = T^ = Ecnn" (5.1.3) 
^" n=0 

As one approaches a generic complex singular point denoted sing from within the complex- 
temperature extension of the FM phase, x is assumed to have the leading singularity 

X ~ ^sing{^ - Ul UsingY^''"'<> {l +ai(l -u/Using) + •••) (5.1.4) 

where A'^^^^ and 7^^^^ denote, respectively, the critical amplitude and the corresponding 
critical exponent, and the ... represent analytic confluent corrections. One may observe that 
we have not included non-analytic confluent corrections to the scaling form in eq. ( |5.1.4| ). 



The reason is that, as discussed in our earlier work |T^, previous studies have indicated that 



they are very weak or absent for the 2D Ising model. We proceed to our results. 
5.2 Singularity at w = We 

Our dlog Fade results relevant to the singularity in x{u,h = m jl) at u = Ue are given 
in Table 2. Because of the length of our differential approximant results relevant to the 
singularity in x at m = Me, we list these in the Appendix. From these results, we infer the 
location of the singularity to be 

Using. = -0.17157 ± 0.00001 (5.2.1) 

This, together with our knowledge of the exact location of the endpoint of the singular 
line segment (|2.33|) , supports the conclusion that the exact location of this singularity is at 
u = Ue = —0.1715729... given in (|2.34|) . Accepting this conclusion, we plot the values of the 
corresponding exponent 7^ for the differential approximants as functions of the normalised 
distance from this point, i.e., \using. — Ue\/\ue\ and extrapolate to zero distance. (This ex- 
trapolation method is similar to the use of biased differential approximants; in both of these 
approaches, one uses one's knowledge of the exact position of the singularity.) From our 
extrapolation, we obtain 

7^ = 1.25 ±0.01 (5.2.2) 

This strongly supports the following inference for the exact value of this exponent, which we 
shall make: 

< = l (5.2.3) 
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r AT- / nl 

[N/D\ 


losing 


1 I/ 1 / U; 1 

1 "'Sing "e|/ | "e| 


7e 


[5/3] 


-0.1716072 


2.0 X 10 * 


1.2492 


[4/4] 


-0.1715705 


1.4 X 10~^ 


1.2464 


[5/4] 


-0.1715636 


5.4 X 10"^ 


1.2459 


[6/4] 


-0.1715646 


4.8 X 10"^ 


1.2460 


[3/5] 


-0.1715802 


4.3 X 10"^ 


1.2471 


[4/5] 


-0.1715628 


5.8 X 10~^ 


1.2458 


[5/5] 


-0.1715615 




1.9 X 10" ' 


1.2160 


[6/5] 


-0.1715598 


7.6 X 10"-'' 


1.2457 


[7/5] 


-0.1715670 


3.4 X 10"^ 


1.2462 


[4/6] 


-0.1715647 


4.8 X 10"^ 


1.2460 


[5/6] 


-0.1715674 


3.2 X 10"^ 


1.2462 


[6/6] 


-0.1715667 


3.6 X 10"^ 


1.2462 


[7/6] 


-0.1715667 


3.6 X 10~^ 


1.2462 


[8/6] 


-0.1715665 


3.7 X 10"^ 


1.2462 


[5/7] 


-0.1715668 


3.6 X 10"^ 


1.2462 


r / — -| 

[6/7] 


-0.1715667 


3.6 X 10~^ 


1.2462 


[7/7] 


-0.1715670 


3.4 X 10"^ 


1.2462 


[8/7] 


-0.1715664 


3.8 X 10~^ 


1.2461 


[9/7] 


-0.1715664 


3.7 X 10~^ 


1.2461 


[6/8] 


-0.1715672 


3.3 X 10"^ 


1.2462 


[7/8] 


-0.1715664 


3.8 X 10"^ 


1.2461 


[8/8] 


-0.1715665 


3.7 X lO"'' 


1.2461 


[9/8] 


-0.1715754 


1.5 X 10~^ 


1.2501 


[10/8] 


-0.1715673 


3.2 X 10~^ 


1.2463 


[7/9] 


-0.1715665 


3.7 X 10"^ 


1.2461 


[8/9] 


-0.1715666 


3.6 X 10"^ 


1.2462 


fQ /Ql 

[9/9J 


— U. i ( ioo ( o 


O.z X iU 


1 0/|(3Q 

i .Z40o 


[10/9] 


-0.1715677 


3.0 X 10-5 


1.2464 


[11/9] 


-0.1715696 


1.9 X 10-^ 


1.2470 


[8/10] 


-0.1715674 


3.2 X 10-5 


1.2463 


[9/10] 


-0.1715686 


2.5 X 10-5 


1.2466 


[10/10] 


-0.1715699 


1.7 X 10-5 


1.2471 


[9/11] 


-0.1715701 


1.6 X 10-5 


1.2472 



Table 2: Values of pole near = —(3 — 2^/^) = —0.171572875..., normalised distance from 

this point, \using—Ue\/\ue\, and exponent 7^ from dlog Fade approximants to low-temperature 

series for Xr{u, h = i7i/2), starting with the series to 0{u^). We list only the approximants 

which satisfy the accuracy criterion \using — Ue\/\ue\ < 2 x 10"^. 
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I AT / Til 

[N/D\ 


"Using 


1 -1 1 

1 Using -L | 


7i 


[4/5] 


-1 r\ r\ f~\ A A A A 

1.0024444 


2.4 X 10"^ 


2.3646 


[6/6] 


0.9894269 


1.1 X 10 


2.3347 


[5/7] 


0.9920029 


8.0 X 10"^ 


2.3623 


[6/7] 


0.9817550 


1.8 X 10"^ 


2.2642 


r — / — n 

[7/7] 


0.9878026 


1.2 X 10"^ 


2.3213 


[8/7] 


0.9988081 


1.2 X 10"^ 


2.4483 


[6/8] 


0.9892965 


1.1 X 10"^ 


2.3380 


[8/8] 


0.9922151 


7.8 X 10"^ 


2.3641 


[9/8] 


0.9969711 


3.0 X 10"^ 


2.4255 


[10/8] 


0.9974836 


2.5 X 10"^ 


2.4332 


1 r\'\ 

[7/9] 


0.9934910 


6.5 X 10"^ 


2.3814 


Pr-i 1 r\'\ 

[8/9] 


1.0054945 


5.5 X 10"^ 


2.5876 


f9/9l 


9975390 


2.5 X 10"^ 


2.4342 


[10/9] 


0.9967646 


3.2 X 10"^ 


2.4231 


[11/9] 


0.9982017 


1.8 X 10-3 


2.4446 


[8/10] 


0.9980098 


2.0 X 10-3 


2.4427 


[9/10] 


0.9992196 


7.8 X 10-^ 


2.4648 


[10/10] 


0.9984417 


1.6 X 10-3 


2.4495 


[9/11] 


0.9985604 


1.4 X 10-3 


2.4519 



Table 3: Values of pole near m = 1, normalised distance from this point, \using — 1|, and 
exponent from [M/N] dlog Pade approximants to low-temperature series for Xr- We list 
only the differential approximants which satisfy the accuracy criterion \using — 1| < 2 x 10^^. 

To calculate the critical amplitude for x at the u = Ue singularity, we use the standard 
method of analysing Pade approximants to the series {—XrY^'^'' (where the minus sign is 
inserted because cq = — 1 in eq. (|5.1.2| )). For A'^ as defined in eq. ( ^.1.4|) with Using = Me, 
we obtain 

A'^ = -0.11515 ± 0.00020 (5.2.4) 
5.3 Singularity at w = 1 

Our dlog Pade results for Xr relevant to the singularity at u = 1 are given in Table 3. We 
obtain more precise results from differential approximants; these are given in Table 4 and 
its continuation. From these results, we obtain the location of the singularity as 
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[L/Mo;Mi] 




1 '^sing 1 1 


7i 


[0/7; 7] 


1.0019932 


2.0 X 10-3 


2.5262 


[0/7; 9] 


0.9983638 


1.6 X 10-=^ 


2.4478 


[0/8; 9] 


0.9974253 


2.6 X 10-3 


2.4337 


[0/8; 10] 


0.9981341 


1.9 X 10^3 


2.4447 


[0/9; 9] 


0.9999105 


0.89 X 10^^ 


2.4776 


[0/9; 10] 


0.9984608 


1.5 X 10-3 


2.4502 


[0/10; 9] 


0.9982548 


1.7 X 10-3 


2.4461 


[1/4 


5] 


1.0013165 


1.3 X 10-3 


2.4393 


[1/6 


8] 


0.9989524 


1.0 X 10-3 


2.4554 


[1/7 


9] 


0.9981225 


1.9 X 10-3 


2.4445 


[1/8 


8] 


1.0019138 


1.9 X 10-3 


2.5153 


[1/8 


9] 


0.9982128 


1.8 X 10-3 


2.4460 


[1/8; 


10] 


0.9981111 


1.9 X 10-3 


2.4444 


[1/9 


9] 


0.9990590 


0.94 X 10-3 


2.4617 


[2/4 


6] 


1.0004562 


0.46 X 10-3 


2.4545 


[2/6 


4] 


1.0008071 


0.81 X 10^3 


2.4534 


[2/6 


6] 


1.0012170 


1.2 X 10-3 


2.4958 


[2/6 


7] 


0.9997299 


2.7 X 10-^ 


2.4678 


[2/6 


8] 


0.9973495 


2.7 X 10-3 


2.4543 


[2/7 


6] 


0.9996006 


4.0 X 10-^ 


2.4651 


[2/7 


9] 


0.9981442 


1.9 X 10-3 


2.4450 


[2/8 


8] 


1.0006495 


0.65 X 10-3 


2.4969 


[2/8 


9] 


0.9985512 


1.4 X 10-3 


2.4513 


[2/9 


8] 


0.9983746 


1.6 X 10-3 


2.4474 


[3/6 


6] 


1.0005153 


0.52 X 10-3 


2.4841 


[3/6 


7] 


0.9973962 


2.6 X 10-3 


2.4232 


[3/6 


8] 


0.9974009 


2.6 X 10-3 


2.4201 


[3/7 


6] 


0.9971786 


2.8 X 10-3 


2.4189 


[3/7 


7] 


0.9977559 


2.2 X 10-3 


2.4302 


[3/7 


8] 


0.9981732 


1.8 X 10-3 


2.4407 


[3/7 


9] 


0.9980242 


2.0 X 10-3 


2.4651 


[3/8 


6] 


0.9975202 


2.5 X 10-3 


2. 1227 



Table 4: Values of pole near u — 1, normalised distance from this point, \using — M, and 
exponent 'y[ from differential approximants to low-temperature series for Xr- We list only 
the differential approximants which satisfy the accuracy criterion \using — 1| < 2 x 10-3. 
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[3/8 


7] 


0.9980148 


2.0 X 10-^ 


2.4375 


[3/8 


8] 


0.9986544 


1.3 X 10-3 


2.4508 


[3/9 


7] 


0.9991902 


0.81 X 10-3 


2.4621 


[4/4 


4] 


1.0024665 


2.5 X 10-3 


2.4681 


[4/5 


7] 


0.9979715 


2.0 X 10-3 


2.4379 


[4/6 


8] 


0.9998460 


1.5 X 10-^ 


2.4832 


[4/7 


5] 


0.9982928 


1.7 X 10-3 


2.1505 


[4/7 


7] 


0.9989074 


1.1 X 10-3 


2.4584 


[4/7 


8] 


0.9990659 


0.93 X 10-3 


2.4624 


[4/8 


7] 


0.9990450 


0.96 X 10-3 


2.4618 


[5/4 


5] 


1.0009762 


0.98 X 10-3 


2.4429 


[5/5 


7] 


0.9987068 


1.3 X 10-3 


2.4543 


[5/6 


7] 


0.9993881 


0.61 X 10-3 


2.4709 


[5/6 


8] 


0.9991317 


0.87 X 10-3 


2.4644 


[5/7 


6] 


1.0014642 


1.5 X 10-3 


2.5280 


[5/7 


7] 


0.9990811 


0.92 X 10-3 


2.4629 


[5/8 


6] 


1.0000180 


1.8 X 10-^ 


2.4888 


[6/4 


6] 


0.9993653 


0.63 X 10-3 


2.5040 


[6/5 


7] 


1.0002641 


2.6 X lO-'' 


2.4944 


[6/6 


7] 


0.9991189 


0.88 X 10-3 


2.4641 


[6/7 


5] 


1.0003178 


3.2 X 10-^ 


2.4956 


[6/7 


6] 


0.9993080 


0.69 X 10-3 


2.4698 


[7/5 


6] 


0.9990014 


1.0 X 10-3 


2.4565 


[7/6 


5] 


0.9983233 


1.7 X 10-3 


2.4405 


[8/4 


6] 


0.9993253 


0.67 X 10-3 


2.4757 


[8/5 


5] 


0.9987939 


1.2 X 10-3 


2.4448 


[8/5 


6] 


0.9994784 


0.52 X 10-3 


2.4852 


[8/6 


5] 


0.9981646 


1.8 X 10-3 


2.4390 


[9/4 


6] 


0.9986630 


1.3 X 10-3 


2.4600 



Table 5: Continuation of table of differential approximants for the singularity in at m = 1. 
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= 0.999 ± 0.001 (5.3.1) 

From this and our determination of the phase boundaries ( p.32D -( p.33| ), we infer that the 
exact location of this singularity is at m = 1. Given this conclusion, we then plot the values 
of 7^ as a function of the distance from u = 1 and extrapolate to zero distance. This yields 
the value 

7i = 2.50 ±0.01 (5.3.2) 

This strongly supports the following inference that we shall make for the exact value of this 
exponent: 

< = l (5.3.3) 

Hence, in particular, 

7; = 27^ (5.3.4) 

We note that the relation ( |5.3.4| ) can be understood if one re-expresses x as a function of the 
elliptic modulus variable /c< in eq. ( |2.7| ), since /c< diverges at m = 1 with an exponent which 
is twice as large as the exponent describing its divergence at m = Ue, i.e., /c< ~ (1 — u)~^ as 
u 1, while fc< ~ (1 — u/wg)"^/^ as u ^ u^- 



5.4 Singularity at u = —1 

We have studied the singularity in x at m = — 1 (denoted Us) by using the series for x in the 
variable u and also transforming this series to one in the elliptic modulus variable /c<. The 
series in A;< showed a greater sensitivity to this singularity, and therefore we concentrate on 
the results from our analysis of this series. The reason for this greater sensitivity is clear; 
the series in u is strongly affected by the fact that, as one can see from Fig. 2(a), there is 
an intervening singular line segment protruding into the FM phase and ending at m = Wg, 
in front of the point u = — 1 as one moves out from the origin along the negative Re{u) 
axis. The transformation from u to A;< maps the singular endpoint at m = Me away to —00 
and maps the singularity at u = 1 to 00 in the fc< plane, thereby leaving the singularity at 
M = — 1 as the nearest to the origin. Specifically, the image of the singular line segment from 
u = Ue leftward to m = —1 is the semi-infinite line segment from —00 to —1 in the A;< plane. 
The line segment from u = 1/ue to u = —1 has the same image, again the segment from 
—00 to —1 in the /c< plane, while the unit circle |m| = 1 maps to the line segment from 1 to 
00 in this plane. The series in fc< has the form x = (l/4)(/c<)^ I]^o '^n(^<)"' ^^^^ before, 
we actually analyse the reduced function Xr = 4(A;<)^^x. Using the Taylor series expansion 
of /?< near u = —1, 

fc< = -1 - 2-5(1 + u)^ + 0((1 + uf) (5.4.1) 
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it follows that as fc< —1 and u —>■ —1, the singular form x ~ (1 + corresponds 
to X ~ (1 + with 

Our results from the differential approximants for this series are given in Table 6 and its 
continuation. Because the actual pole positions have small imaginary parts, typically a few 
times 10~^ of the size of the real part, there are resultant imaginary parts in the values of the 
corresponding exponent 7^ from the differential approximants. Since the exact singularity 
in xi"^) is at the real value u = —1, and since h = m/2) is real for real m, we know that 
7^ at u = (A;<) = — 1 is real. Given this and the relation ( |5.4.2| ), it follows that we may take 
only the real parts of the exponents from the differential approximants to the series in 
and we do so. 

From this study, we obtain for the position of the singularity 

{k<)sing = -0.9993 ± 0.0001 (5.4.3) 

consistent with the expectation {k^)sing = —1, or equivalently. Using = Ug = —1. As is 
evident, the values of i?e(7^) are almost all slightly below 0.25; however, when we carry out 
our method of plotting the values as a function of the distance |(fc<) + 1| and extrapolating 
to zero distance from the exact singularity, the extrapolated value is actually slightly above 
0.25. Accordingly, we give a conservative estimate 

7:,,< = 0.250 ± 0.020 (5.4.4) 

and hence, using ( p.4.2|) , 

7; = 1.00 ±0.08 (5.4.5) 
This supports the conclusion, which we shall draw, that the exact value of this exponent is 

7: = 1 (5.4.6) 

We show our summary of exponents in table 8. The exponent relation C('^^ph. + '^l3'u~^l'u,ph. = 2 
is evidently satisfied at all three of the singularities u = Ue, u = 1, and u = —1. 

6 Extraction and Analysis of Low- Temperature Series 

for x^^^ 

We have also investigated the complex-temperature singularities in the staggered suscep- 
tibility x''"^ for the present model. To do this, we have extracted and analysed the low- 
temperature series expansion for this function, using the low-temperature, high staggered 
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7.3 X 10-4 
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-0.9992847 + 0.0002823i 
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Table 6: Values of pole {k<)sing near A;< = —1, i.e. u = —1 (denoted Ug), distance from this 
point, \{k^)sing + 1|, and real part of exponent, Re{'~fgf^^) from differential approximants to 
low-temperature series for Xr(k^,h = m/2). We list only the approximants which satisfy 
the accuracy criterion \{k^)sing + 1| < 1 x 10-^. 
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0.2218 
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-0.9995104 + 0.000392H 


6.3 X 10"^ 


0.2390 


[10/4; 4] 


-0.9994550 - 0.0001668i 


5.7 X 10"^ 


0.2311 
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-0.9991991 + 0.0001605i 
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0.2164 


[10/5; 4] 


-0.9992330 + 0.0001856i 
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Table 7: Continuation of table of differential approximants for Xr series in the variable k. 
near u — k< — —1. 



u 


^u,FM 


^u,AFM 






^u,ph. + 2/3^ + j'u^pfi. 


lie = -(3 -2^/^) 


1 




-1/8 


5/4 


2 


1 


finite* 


finite* 


-1/4 


5/2 


2 


Us = -1 


finite 


finite 


1/2 


1 


2 



Tabic 8: Exponents at singularities in the 2D Ising model with h = ±i7r/2. The results for 
a'^ and Pu are exact; the results for 7^ are our conclusions for the exact values from our series 
analysis. The notation — indicates that the point cannot be approached from within the 
given phase. For the singularity of C at w = 1 (marked with a *), the values of a' correspond 
to evaluating K = — (l/4)ln(l) = on the principal Reimann sheet of the logarithm, as 
discussed in the text. 
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field series expansions for tfie free energy of the Ising model on the square lattice calculated 



by the King's College group These are denoted antiferromagnetic polynomials in 



these papers and were calculated to order 0{w^^) in Ref. pSf, where, as before, w = 1/u is 



the low-temperature expansion variable in the AFM phase. The antiferromagnetic polyno- 
mials were apparently not calculated to higher order subsequently We have extracted 



from these the resultant low-temperature series expansion for x^""^ for h = in/2, which is 
X^'') {h = m/2) = 4w;2 [ - 1 - 8«;2 + 2Aw^ - 135«;^ + QA&w^ - 3336w^ 

+17240w^ - 90501^;^ + 479192^^ + 0{w^^)\ (6.1) 

For reference, we recall that the series for x^*^) for /i = on this lattice is 13^, |25[ 



X^"^)(h = 0) = 4w=^[i+4w2+8w^+39m;^+152w^+672w^+3016w^+13989w^+66664w9+O(w1°) 

(6.2)" 

The series ( |6.1[ ) is much shorter than the one which we extracted for x, given by eq. ( |5.1.2[ ) 
and Table 1, and hence one does not expect to derive results for x*^"^ which are as precise 
as those which we obtained for x- As before, we have used both dlog Fade and differential 
approximants for this analysis. 

We study first the vicinity of the singular point w = 1. In Table 9 we list the pole locations 
and resultant exponents for the approximants which satisfy the accuracy requirement \wsing — 
1| < 0.05. One would not normally expect a dlog Fade approximant of such low order as [1/2] 
to yield an accurate result; however, it happens that the denominator of this approximant 
is oc (1 — w){l — (ll/2)i(7), so that it produces a location for the pole which is exact. 
For this reason, it yields a much better determination of the associated exponent than 
would otherwise have been the case. This fortuitously accurate approximant, and the best 
differential approximant, [3/2; 2], both give 7^^^ values of about 2.5. From the full set, we 
infer the crude result 

7i,, = 2.5±0.5 (6.3) 

This is consistent with the exact value '~^[^ = 5/2 and hence with the equality 7^^ ^ = -y'l- 
However, clearly the results for 'y[ ^ are much less precise than our determination of 7J. 

We also studied the series in the vicinity of the singular endpoint w = We = —{3 — 2'^/^) 
(i.e., u = Uoe = l/we = — (3 + 2^/^)). To optimise the sensitivity, we calculated and analysed 
series in transformed variables to map the singularity at w = 1 away. We required these 
variables to be equal to w for small w and to map w = ±00 to ±00, respectively. Two 
such variables were w' = w{l + w/8){l — w)~^ and w" = w{l — w)^^smh.w. The series 
in the transformed variables did slightly better in locating the pole positions in w'^ and w" 
corresponding to We- The dlog Fade and differential approximants indicated that x^"^ has a 
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[0/3; 4] 


0.955169 


4.5 X 10-2 


2.057 


[0/4; 2] 


0.955792 


4.4 X 10-2 


2.172 


[2/2; 2] 


0.999046 


0.95 X 10-^ 


2.554 


[3/2; 2] 


0.952282 
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Table 9: Values of pole near w = 1 (denoted Wi), normalised distance from this point, \wsing — 
1|, and exponent 'j[ „ from dlog Pade and differential approximants to low-temperature series 
for xi'^\ We list only the Pade and differential approximants which satisfy the accuracy 
criterion \wsina — 1| < 0.05. 



divergent singularity at We and yielded values for the associated exponent y^^ g^ in the range 
from about 0.2 to 0.4. Given our exact results a^g = 1 for the specific heat and Poe = — 1/4 
for the staggered magnetisation, a value within the above range for ^ would indicate a 
violation of the exponent relation a'g^ + 2Poe + 'loe = 2. In this context, it is of interest to note 
that we have already found violations of the relation a + 2/5 + 7 = 2 at complex-temperature 
singularities, e.g., in the zero-field Ising model on the square lattice at u = Ug = —1, as 
approached from within the PM phase, where = 0, /^^ = 1/4, and 7^ < (since x has a 



finite non-analyticity for the approach from within the PM phase) [18|, and in the zero-field 



Ising model on the honeycomb lattice, at the point z = Zf = —1, a.s approached from within 



the FM phase, where a'f, = 2, f3e = -1/4, and 7^ = 5/2, so that + 2/?^ + 7^ = 4 pO 



7 Complex- Temperature Behaviour of the Correlation 
Length 

In this section we shall study the complex-temperature behaviour of the correlation length. 
To do this, we make use of a calculation of the asymptotic form of the spin-spin correlation 
function along a row (or equivalently, column), < cro,oC"n,o for large n [Q (where, without 
loss of generality, one may take n > 0. From this calculation, carrying out an analytic 
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continuation to complex temperature, we obtain, for n oo, 

< OQfiOnfi >conn. {2 / 'r){1 - U^Y^ n-^u{-uY 

= -(2/7r)(l - u)-^'\l + 6m + u'Y^/^n'^ui-uY (7.1) 

where we have used the exact expression for M, (f4.1|) . This analytic continuation ap- 
plies within the FM phase. Extracting the correlation length ^ in the usual way as = 
— lim^^oo '''"^ ln(< cToCTr >conn.), whcrc r = |r|, we find 

CoL = -M-«) (7.2) 

For usual physical second-order critical points, one can use the connected 2-spin correlation 
function for any r, with |r| oo, to extract the correlation length ^. However, in our 
previous work |18|, |l9l, we found that at the complex-temperature singular point u = Us = —1 
in the zero-field Ising model on the square lattice, the correlation length defined from the 
diagonal connected 2-spin correlation function diverges with a different exponent, I's^diag = 2, 
than the exponent z/^ = 1 describing the divergence in the correlation length defined from 
off-diagonal (e.g., row) correlation functions. In view of this, we include the suffix row in 



( [r.2|) for clarity. We now consider three particular singular points which can be approached 



from within the complex-temperature FM phase, viz., u = Ue, u = —1, and u = 1. As 
u Ue, the 2-spin correlation function ( [7111) diverges, like (1 — u/ue)~^^'^, because of the 
divergence in the prefactor M^, but the correlation length ^row remains finite, with = 
— ln(— Me) = 1.7627... at M = Ue- If this feature of a finite correlation length applied to 
all of the connected 2-spin correlation functions, precisely at u = Ue as well as for points 
approaching Ue from within the complex-temperature FM phase, then by the same argument 
as was used in Ref. |^, it would follow that the only singularity in x would arise from the 
divergent prefactor.0 We know, however, that the above premise cannot be true, since 
then the susceptibility exponent at Ue (as approached from within the FM phase) would 
be 1/4, whereas we found that j'e — 5/4. The fact that the susceptibility diverges with an 
exponent different from that arising from the divergent prefactor shows that at least 
some connected 2-spin correlation functions must decay like a power law, i.e. the associated 
correlation length must be divergent, at m = Ue- To get more information on this, it would 
be useful to carry out analytic calculations of the asymptotic forms of the general 2-spin 
correlation functions < crooCmn > in the present model, near to and at this singular point. 



■^Define < (Tofr > = M^c(r). Then x = J2r '-('")■ purposes of analysing divergences, the asymp- 
totic behaviour of the sum can be approximated by that of the integral J (Pr c(r) for large r (with a 
short-distance cutoff on the latter). If c(r) ~ r~Pe~^^^ as r ^ oo, this integral is finite. Therefore, a 
divergence in x would arise solely from a divergence in the prefactor M^. This was noted in Ref. []l9|. 
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As M 1, the 2-spin correlation function (|7. 1| ) again diverges, like (1 — m)"^/^, and 
the correlation length ^row is finite: = — ln(— 1) = —in (for the principal Riemann 

sheet of the logarithm). This is a case similar that discussed in |^ where i?e(^~^) = but 
Jm(r') 7^ 0. 

As M ^ — 1, each 2-spin correlation function is finite, but the correlation length does 
diverge, with exponent 

K = 1 (7.3) 

If one were to use the exponent relation 7^ = z/g(2 — 77^), then with our inference 7^ = 1 in eq. 
( p.4.6| ), it would follow that 77^ = 1. However, we have shown previously |]18| that one must 



use caution in trying to apply such exponent relations at complex-temperature singularities, 
since different connected spin-spin correlation functions may be characterised by correlation 
lengths which diverge with different exponents u. 

This type of analysis can also be done with the staggered 2-spin correlation functions 

(c.f. (EnD) 

< 5"o,oO"ra,0 >conn. = (^1)" < O"0,0O"n,0 >conn. (7-4) 

In particular, as m — >■ 1/ue (i.e., w Wg), these correlation functions diverge, like (1 — 
w/we)~^^^, because of the divergent prefactor M^j. However, the correlation length ^row,AFM 
remains finite. If this behaviour characterised all of the staggered 2-spin correlation functions, 
at We as well as in the vicinity of We, then the only divergence in x^"^ would arise from the 
prefactor, and hence 7^^ = 1/4. This value is consistent with our results from the 
analysis of the low-temperature series for x^"'\ This merits further study. 

8 Exact Solution at u = 1 for Arbitrary H 

In the body of this paper, we have investigated singularities in the square lattice Ising model 
as functions of complex-temperature, for the fixed value of external magnetic field, h = m/2 
(or h = —m/2). It is also of interest to study the complementary problem of singularities 
as a function of h for fixed K or u. Indeed, in pursuing such a study, Yang and Lee were 
led to their celebrated circle theorem on the zeros of the partition function for the Ising 
model in the complex e^^ plane 32|. Here, we would like to mention some elementary 



results which elucidate how various quantities become singular at a particularly simple point, 
M = 1, as /i is varied. These results may be combined with our determination of the exact 
singularities in /, t/, C, and M as one approaches this point by varying u. At J = 0, hence 
K = Q and u = 1, the partition function reduces to a single-site problem, which can easily be 
calculated exactly for arbitrary if, dimensionality, and lattice type. We find, independent 
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of the dimensionality and lattice type, 



f{u = l,h) = ln(2cosh/i) 
U = -Htanhh 



^ cosh^ h 
M = tanh h 
1 

X 



8.1) 

8.2) 

8.3) 
8.4) 

cosh^ h ^ ^ 

Further, the m-point correlation functions factorise trivially and are independent of the 
positions of the spins: 

< ^ri ■ • ■ fXr^ > = < (Jr, = M'" = [tauh (8.6) 

To study the singularities of these functions, we must define new exponents, since the usual 
critical exponents apply to singularities of thermodynamic quantities as functions of T. To 
avoid a profusion of new symbols, we shall use the same Greek letters as for the respective T- 
dependent singularities in thermodynamic quantities, but use a superscript (h) to indicate 
that they describe the singularity as a function of h for K = 0. Thus, for the leading 
singularity in the specific heat, as a function of h, at the point h = hs {s denotes a generic 
singularity here) for fixed K = (hence u = 1), we shall write 



C(/.)..„,.~4"U.(/^-/^.)""''^- (8.7) 

where dir. denotes the direction, in the complex h plane, from which one approaches the 
singular point hg. Similarly, we shall write 



(h) 



V • ~ A^^^ (h- h V^^d^r. 

Xsmg. ^x,s,dir.\"' "'sj 



(h) 



and so forth for the singularities in other quantities. From (|8.3|) , it is evident that C diverges 
for 

h = {2n + l)— , n^Z (8.9) 

with corresponding exponent 

af^ = 2 (8.10) 
for any direction of approach to any of the singular points (|8.9|). The internal energy itself 



also diverges at these points, with the exponent — 1 = 1 and vanishes at the set 

h = nin , n E Z (8-11) 
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The magnetisation vanishes and diverges at the same set of points as the internal energy U 



(c.f. eqs. and ( p.9|) with the respective exponents 

(^i%ro = 1 (8-12) 

and 

& = -1 (8.13) 

again, independent of the direction of approach to these points in the complex h plane. The 
susceptibility diverges at the points with exponent 

7f ) = 2 (8.14) 

It is useful to evaluate the general results above for the interesting special case of complex 
h = hr + m /2 where hr is real. For this case, we have 

/ = ln(2i sinh /i,.) (8.15) 

U = ^ (8.16) 

tanh hr 



smh n. 



M = -^TT (8-18) 
tanh hr 

X = -- (8-19) 
smh hr 

and 

< f^ri • ■ ■ f^r^ > = (tanh/i^,)-'" (8.20) 

As /ir — > SO that h ^ 111/2, we recover our previous results, that U and M diverge linearly 
while C and x diverge quadratically. 

For loose-packed lattices, it is straightforward to extend these results to consider a stag- 
gered rather than uniform external field, Hgt, i.e. to consider the partition function 

Z=Y^ exp(5:(-l)P(")/i,ia„) (8.21) 

where p{n) was defined in eq. ( [4.12 ) and h^t = jSHgt- Since the summations over the spins 
on the even and odd sublattices are decoupled, one finds the same equations as before, but 
with H replaced by Hgt, i.e., / = ln(2 cosh /i^j), U = — Hgt tanh. hgt, C = kshl^/ cosh^/i^t, 
and Mst = ianhhst- The staggered susceptibility is x''"'' = 1/ cosh^/i<jt. 



31 



9 Conclusions 



In this paper, we have studied a natural generahsation of an exactly solved model from real 
non-negative temperature to complex temperature. This is the Ising model on the square 
lattice in an external magnetic field given by (3H = ±i7r/2, first solved by Lee and Yang 
We have worked out the complex-temperature phase boundaries, as shown in Fig. 2. 
We have also extracted a low-temperature series expansion for the susceptibility, x- From 
an analysis of this series using dlog Fade and differential approximants, we conclude that 
X has divergent singularities at m = Me = — (3 — 2'^/^) with exponent 7^ = 5/4, at -u = 1 
with exponent 'j[ = 5/2, and at m = = — 1 with exponent 7^ = 1. We have also 
studied the staggered susceptibility. Using exact results, we have determined the complex- 
temperature singularities of the specific heat and the uniform and staggered magnetisation. 
We are currently in the process of extending our studies to other 2D lattices. The findings 
show again that even though the Ising model has a very simple Hamiltonian, it exhibits a 
fascinating richness of properties. 

This research was supported in part by the NSF grant FHY-93-09888. One of us (R.S.) 
thanks Frof. C. N. Yang for a discussion of Ref. and Frofs. D. S. Gaunt and A. J. 
Guttmann for discussions of the current status of low-temperature series expansions. 
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Appendix 

Here we list our differential approximant results relevant for the singularity in h ~ 

m/l) at M = Ue- 
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[VA^o;Mi] 




1 1 /I 1 

1 "'sing "e|/ | "e| 



7e 


[0/4 


6] 


-0.1715614 


6.7 X 10~^ 


1.2455 


[0/5 


5] 


-0.1715743 


8.1 X lO"*" 


1.2474 


[0/5 


6] 
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7.2 X 10"^ 


1.2470 


[0/5 


7] 
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1.9 X 10"^ 
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[0/6 


4] 


-0.1715688 


2.4 X 10~^ 


1.2466 


[0/6 


5] 


-0.1715717 


6.6 X 10~^ 


1.2470 


[0/6 


7] 
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[0/6 


8] 
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4.1 X 10"^ 
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[0/7 


5] 
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3.9 X 10"^ 
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7] 


-0.1715659 


4.1 X 10"^ 


1.2460 


[0/7 


8] 


-0.1715660 


4.0 X 10~^ 


1.2460 


[0/7 


9] 


-0.1715628 


5.8 X 10"^ 


1.2456 


[0/8 


6] 


-0.1715658 


4.1 X 10"^ 


1.2460 


[0/8 


7] 


-0.1715660 


4.0 X 10"^ 
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8] 


-0.1715654 


4.3 X 10"^ 
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[0/8; 10] 


-0.1715690 


2.2 X 10"^ 


1.2469 


[0/9; 7] 


-0.1715599 


7.6 X 10~^ 


1.2456 


[0/9; 8] 


-0.1715549 


1.0 X 10"^ 


1.2468 


[0/9; 9] 


-0.1715754 


1.5 X 10"^ 
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[0/9; 10] 


-0.1715716 


7.2 X 10^*^ 
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[0/10; 8] 


-0.1715584 


8.5 X 10"^ 


1.2456 


[0/10; 9] 
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[1/5 


5] 


-0.1715621 


6.3 X 10~^ 


1.2449 


[1/5 


6] 


-0.1715609 


7.0 X 10"^ 


1.2447 


[1/5 


7] 


— U. i ( iDDZD 


O.U X iU 




[1/6 


5] 


-0.1715609 


7.0 X 10-5 


1.2446 


[1/6 


6] 


-0.1715615 


6.7 X 10-5 


1.2448 


[1/6 


7] 


-0.1715649 


4.7 X 10-^ 


1.2457 


[1/6 


8] 


-0.1715676 


3.1 X 10-5 


1.2466 


[1/7 


5] 


-0.1715624 


6.1 X 10-5 


1.2450 


[1/7 


6] 


-0.1715647 


4.8 X 10-5 


1.2457 



Table 10: Values of pole near Ue = -(3 - 2^/2) = -0.171572876..., normalised distance 
from this point, \using — Ue\/\ue\, and exponent 7g from differential approximants (DA's) to 
low-temperature series for Xr- We list only the differential approximants which satisfy the 
accuracy criterion \using — Ue\/\ue\ < 1 x 10~^. 
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[1/7 


i—r\ 
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-0.1715649 
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[2/7 


6] 


-0.1715657 
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_r; 
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1.2467 
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o 
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[2/7 
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[2/8 


8] 


— U.i / iODoD 


Z.O X iU 


-5 
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[2/8 
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[2/9 
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-0.1715693 


2.1 X 10' 


-5 


1.2468 


[2/9 


8] 


-0.1715785 


3.3 X 10' 


~5 


1.2506 


[3/4 


4] 


-0.1715556 


1.0 X 10' 


-4 


1.2434 


[3/4 


5] 


-0.1715619 


6.4 X 10' 


-5 


1.2449 


[3/4 


6] 


-0.1715632 


5.6 X 10' 


-5 


1.2452 



Table 11: Continuation of DA Table for Ug Singularity. 
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[3/5 
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-0.1715674 


3.2 X 10-5 


1.2441 



Table 12: Continuation of DA Table for Ug Singularity. 
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Table 13: Continuation of DA Table for Ug Singularity. 
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[7/5 
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-0.1715652 


4.5 X 10- 


-5 


1.2388 


[7/7 


5] 


-0.1715689 


2.3 X 10- 


-5 


1.2458 


[8/4 


2] 


-0.1715791 


3.6 X 10- 


-5 


1.2505 


[8/4 


3] 


-0.1715670 


3.4 X 10- 


-5 


1.2464 


[8/4 


4] 


-0.1715679 


2.9 X 10- 


-5 


1.2466 


[8/4 


5] 


-0.1715683 


2.7 X 10" 


-5 


1.2465 


[8/4 


6] 


-0.1715609 


7.0 X 10- 


-5 


1.2283 


[8/5 


3] 


-0.1715677 


3.0 X 10" 


-5 


1.2467 


[8/5 


4] 


-0.1715682 


2.7 X 10- 


-5 


1.2465 


[8/5 


5] 


-0.1715691 


2.2 X IQ- 


-r, 


1.2161 


[8/6 


4] 


-0.1715700 


1.7 X 10- 


-5 


1.2476 


[8/6 


5] 


-0.1715621 


6.3 X 10- 


-5 


1.2310 


[9/4 


2] 


-0.1715693 


2.1 X 10- 


-5 


1.2472 


[9/4 


3] 


-0.1715676 


3.1 X 10" 


-5 


1.2467 


[9/4 


4] 


-0.1715684 


2.6 X 10- 


-5 


1.2466 


[9/4 


5] 


-0.1715692 


2.2 X 10- 


-5 


1.2471 


[9/4 


6] 


-0.1715572 


9.1 X 10- 


-5 


1.2152 


[9/5 


3] 


-0.1715658 


4.1 X 10- 


-5 


1.2458 


[9/5 


4] 


-0.1715691 


2.2 X 10- 


-5 


1.2469 


[9/5 


5] 


-0.1715690 


2.3 X 10- 


-5 


1.2460 


[9/6; 4] 


-0.1715709 


1.1 X 10- 


-5 


1.2483 


[10/4 


;2] 


-0.1715630 


5.7 X 10- 


-5 


1.2448 


[10/4 


;3] 
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-5 


1.2473 
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;4] 


-0.1715669 


3.5 X 10- 


-S 


1.2433 


[11/5 


;3] 


-0.1715707 


1.2 X 10" 


-5 


1.2481 



Table 14: Continuation of DA Table for Ug Singularity. 
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Figure Caption 



Fig. 1. Assignment of couplings Ki, i = 1,2,3,4 to links (bonds) of tlie clieckerboard 
lattice. 

Fig. 2. (a) Phases and associated boundaries in the complex u plane for the Ising model 
on the square lattice with h = ±.m/2. The boundaries are given by eqs. (|2.32|) and ( p.33|) 
in the text. In particular, the line segment extends from Ue as given in eq. (|2.34| ) off the 
figure to the left, ending at 1/ue = —(3 + 2^/2) ~ —5.828. FM and AFM refer to phases in 
which M 7^ 0, Mst = and M = 0, Mgt 7^ 0, respectively, (b) Complex-temperature phase 
diagram in the v plane. 

Fig. 3. Path of /t^ corresponding to u = pe*^ for 6* = 7r/6, as p varies from to 30. 
Horizontal and vertical axes are the Re{n'^) and Im{n'^) axes. The image of the singular 
curve ( |2.32|) and line segment ( p.33|) is the dark line from = 1 to = cxd. 
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